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Abstract 



We show that in a large class of supersymmetric models with sponta- 
neously broken B — L symmetry, neutron-antineutron oscillations occur at 
an observable level even though the scale oi B — L breaking is very high, 
vb-l ~ 2 X 10^^ GeV, as suggested by gauge coupling unification and neu- 
trino masses. We illustrate this phenomenon in the context of a recently 
proposed class of seesaw models that solves the strong CP problem and the 
SUSY phase problem using parity symmetry. We obtain an upper limit on 
N — N oscillation time in these models, t^_7v ^ 10^ ~ 10"'^'^ sec. This suggests 
that a modest improvement in the current limit on Tj^_f^ of 0.86 x 10^ sec will 
either lead to the discovery of N — N oscillations, or will considerably restrict 
the allowed parameter space of an interesting class of neutrino mass models. 



I. INTRODUCTION 



It is widely believed that the most natural and appealing explanation of the recent 
neutrino oscillation results is provided by the seesaw mechanism incorporated into ex- 
tensions of the Standard Model that include a local B — L symmetry. The simplest models 
with local B — L symmetry are the left-right symmetric models based on the gauge 
group SU{Z)c X SU{2)l x SU{2)ji x U{1)b-l- These models have the additional virtue 
that they explain the origin of parity violation in weak interactions as a consequence of 
spontaneous symmetry breaking in very much the same way as one explains the strength of 
the weak interaction in the Standard Model. Stability of the Higgs sector under radiative 
corrections calls for weak scale supersymmetry as in the Minimal Supersymmetric Standard 
Model (MSSM). It has recently been shown that if the MSSM is embedded into a left-right 
symmetric framework at a high scale vr ~ 10^^ — 10^^ GeV, as suggested by neutrino os- 
cillation data and by gauge coupling unification, it helps solve some important problems 
faced by the MSSM, viz., the SUSY CP problem [Q], the strong CP problem and the fi 
problem. Supersymmetric models with such a high scale embedding are therefore attractive 
candidates for physics beyond the Standard Model. 

It was noted many years ago that the electric charge formula of the left-right symmet- 
ric models, Q = I3L + hR + allows one to conclude from pure group theoretic arguments 
that parity symmetry breaking implies a breakdown oi B — L symmetry as well with the 
constraint that 2A/3R = —A{B — L). This simple relation is profoundly revealing. It says 
that the neutrinos must be Majorana particles since the lepton number breaking terms in the 
theory must obey |AL| = 2 selection rule. This conclusion follows directly if Higgs triplets 
are used to break SU{2)r symmetry since J3R = 1 for triplets, it also holds when Higgs 
doublets are used for this purpose, since gauge invariance requires the presence of two such 
doublets in the mass term for the neutrinos. Secondly, for purely hadronic baryon number 
violating processes, baryon number must change by at least two units, |Ai?| = 2. This 
means that models based on left-right symmetric gauge structure can lead to the process 
where a neutron transforms itself into an antineutron (A^ — N oscillation 0,1^,0), while they 
may forbid the decay of the proton, which is a AB = 1 process. 

While the above group theory argument predicts the existence of N — N oscillation in 
left-right symmetric models, its strength will depend on the details of the model. Using 
simple dimensional analysis it is easy to find that the lowest dimensional operators that 
contributes to N — N oscillation are six quark operators, a typical one being (u^d'^u^d'^d'^d^) . 
This operator has dimension 9 and therefore the coupling strength scales as Gab=2 ~ jgs, 
where M is the scale of new physics. It is natural to identify M with the scale oi B — L 
(or parity) breaking. The current lower limit on N — N oscillation time, Tj\f_^ > 0.86 X 10« 
sec imphes an upper limit G'ab=2 < 3 x 10~^^ GeV~^. For N — N oscillations to be 
observable then, the scale M should be rather low, M < 10^ GeV. 

One class of models where AB = 2 transition manifests itself through Higgs boson 



This is the direct limit from free neutron oscillation searches. Indirect Hmits which involves some 
reasonable nuclear physics assumptions have been extracted from nucleon decay experiments which 
are slightly more stringent: t^_jv ^ 1-2 x 10^ sec 
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exchange has been discussed in Ref. There it was shown that if the SU{3)c x SU{2)l x 
SU{2)r X U{1)b-l model is embedded into the SU{4:)c x ^^7(2)^ x SU{2)r gauge group, 
then N — N oscillations can arise at an observable level if the SU{4)c breaking scale is in 
the 100 TeV range. In these models, N — N oscillation amplitude is intimately tied to an 
understanding of small neutrino masses via the seesaw mechanism as well as the breaking of 
quark-lepton degeneracy implied by SU{A)c symmetry. The same Higgs field that breaks 
SU{A)c and generates heavy Majorana masses for the right-handed neutrinos also mediate 
N — N oscillations here. With the scale of SU{A)c breaking in the 100 TeV range, these 
models would appear to be incompatible with gauge coupling unification. Furthermore, 
such a low scale of parity breaking would not yield naturally neutrino masses in the range 
suggested by current experiments. If we raise the scale of parity/S't/(4)c breaking to values 
above 10^^ GeV, so that small neutrino masses in the right range are generated naturally, 
then N — N transition amplitude becomes unobservably small in these models.0 

Does the above arguments mean that N — N oscillations are beyond experimental reach 
based on current neutrino oscillation phenomenology? In this Letter we will show that this 
is not the case in a class of attractive seesaw models with local B — L symmetry. We will 
see that in these models a new class of Ai? = 1 operators is induced as a consequence of 
parity breaking. These operators lead to observable N — N oscillation despite the scale vr of 
parity breaking being close to the conventional GUT scale of 2 x 10^^ GeV. In fact, G'ab=2 
increases with vr and therefore one has the inverse phenomenon that increasing vr leads 
to stronger N — N oscillation amplitude. Interestingly, the scale vr implied by neutrino 
masses is such that N — N oscillation should be accessible experimentally with a modest 
improvement in the current limit. We obtain an upper limit of Tj^_j^ < 10^ — 10^*^ sec in 
this class of models. This prediction becomes sharper in a concrete model where flavor 
symmetries reduce considerably the uncertainties in the estimate of Tjv_jv- We emphasize 
that our upper limit is derived in the context of conventional seesaw models of neutrino 
mass without using any special ingredients to enhance N — N oscillation amplitude. This 
should provide new impetus for an improved experimental search for N — N oscillations. 

II. HIGH SCALE SEESAW MODEL AND N - N OSCILLATION 

The basic framework of our model involves the embedding of the MSSM into a minimal 
SUSY left-right gauge structure at a scale vr close to the GUT scale. The electroweak gauge 
group of the model, as already mentioned, is SU{2)l x SU{2)r x U{1)b-l with the standard 
assignment of quarks and leptons - left-handed quarks and leptons (Q, L) transform as dou- 
blets of SU (2) L, while the right-handed conjugate ones (Q'^, L^) are doublets of SU{2)r. The 



counter example where a higher scale of parity violation can go hand in hand with observable 
N - N oscillation was noted in the context of a SUSY SU{2)l x SU{2)r x S'[/(4)c model in Ref. 
|10|. These models possess accidental symmetries that lead to light (~ 100 GeV) diquark Higgs 
bosons even though the scale of parity violation is high. As a result, the N — N oscillation operator 
can have observable strength. Unification of gauge couplings is however difficult to achieve in these 
models. 
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quarks Q transform under the gauge group as (2, 1, 1/3) and Q'^ as (1,2, —1/3), while the 
lepton fields L and L'^ transform as (2, 1, —1) and (1, 2, +1) respectively. The Dirac masses of 
fermions arise through their Yukawa couplings to two Higgs bidoublet $a(2, 2, 0), a = 1 — 2. 
The SU{2)ji X U{1)b-l symmetry is broken down to U{1)y in the supersymmetric limit 
hj B — L = ±.1 doublet scalar fields, the right-handed doublet denoted by x'^(l,2, — 1) ac- 
companied by its left-handed partner x(2, 1, 1). Anomaly cancellation requires the presence 
of their charge conjugate fields as well, denoted as x'^(l,2, 1) and x(2, 1, —1). The vacuum 
expectation values (VEVs) {x^) = {x^) = vr break the left-right symmetry group down to 
the MSSM gauge symmetry. A singlet S is also used to facilitate symmetry breaking in the 
SUSY limit. 

It has recently been shown that if there exists a Z4 R symmetry, the minimal model 
just described will solve the strong CP problem and the SUSY phase problem based on 
parity symmetry 0]. Furthermore, the /i term will have a natural origin. One possible Z4 
assignment was given in Ref. Here we present a slight variant, which yields the same 
superpotential at the renormalizable level as in Ref. and thus preserves all its success. 
Under this Z4, the superpotential W changes sign, as do and The quark fields 
iQiQ^) are even, while {L,L^) transform as The fields {x^XTX^iX'^i^ai S) are all 

odd under Z4. 

The gauge invariant superpotential consistent with this Z4 R symmetry at the renormal- 
izable level is 



This superpotential breaks the gauge symmetry to that of the Standard Model in the SUSY 
limit without leaving any unwanted Goldstone bosons and induces realistic quark masses 
and mxings. 

The baryon number violating processes as well as neutrino masses arise in this model 
from higher dimensional operators induced by Planck scale physics. They will be the main 
focus of the rest of the paper. We shall pay special attention to the relation between the 
neutrino mass and the N — N oscillation time. The relevant dimension four operators in the 
superpotential which are scaled by Mp^ and are allowed by the Z4 symmetry are: 



Operator Oi gives rise to Majorana masses for of order v\/M-p\. Combining this with the 
Dirac neutrino masses arising from Eq. (1), small neutrino masses will be generated by the 
seesaw mechanism. For vr ~ 10^^ — 10^^ GeV, the magnitude of the light neutrino masses 
are in the right range to explain the atmospheric and the solar neutrino oscillation data. 
Operator (^2, which is also invariant under the Z4, leads to baryon number violation. While 
Oi^2 could have their origin in quantum gravity, they may also be induced by integrating 
out vector states that have Z4-inviariant masses of order the Planck scale. 

Note that operators such as L^x'^i LLL'^x'^i QLQ'^x'^ ell's not allowed by the Z4 symmetry. 
If they were present along with O2, they would lead to rapid proton decay. Note also that 
the well known proton decay operator QQQL is not allowed by the Z4 symmetry. In any 



^^S{e'^x'x' + e-'^xx + aS^ - M^) + /^„,Tr($„<l>6)S . 



(1) 



O, = f [{L^xr + (LxY , 

02 = f' [QVVx' + QQQx] 



(2) 
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case its presence would not have been a problem since it is scaled by the Planck mass and 
therefore can lead to a proton lifetime consistent with the present lower limit. 

To see the connection between neutrino masses and the N — N oscillation time Tjv_jv = 
{1/Sm^_f^) qualitatively, first we note that the operator Oi leads to the Majorana mass 

for the right handed neutrino Mr = -p^. The seesaw formula then leads to the relation 

nil, = — T-T^- On the other hand, the operator O2 leads to a AB = 1 operator with strength 
Leaving aside the details of the flavor structure of O2 and how actually 5mj^_fj arises, 
it is clear that we have a simple linear relation between the neutrino masses and the N — N 
oscillation time: 

where C is a dimensional constant which depends only on the details of weak scale physics 
and does not involve the high scale vr. We will evaluate C in the next section. This simple 
relation makes it clear that our present knowledge of the neutrino masses allows a direct 
prediction of the N — N oscillation time in the context of the supersymmetric left-right 
models broken by doublet Higgs fields. 



III. FROM SUPERSYMMETRIC AS = 1 OPERATOR TO N - N OSCILLATIONS 

Let us now proceed to examine the expected N — N oscillation time resulting from 
the AB = 1 operator O2 of Eq. (2). An important point to note here is that since O2 
is a superpotential term with antisymmetric color contraction it must have antisymmetric 
flavor contraction as well. The flavor structure of this operator is then of the type u'^d'^s'^, 
u'^d'^lf or u'^s'^lf in terms of the superfields. We must then use flavor mixings to obtain 
the fermionic operator of the type u^d'^d'^ and then the six quark N — N operator. The 
dominant contribution to this process comes from the Feynman diagram shown in Fig. 1 
which proceeds through the exchange of a gluino and squarks [|TT| and involves two d'^ — b*^ 



mixings. The strength of the AB = 2 operator resulting from Fig. 1 can be estimated to be 



Gab=2 ^ ' , (4) 



where Mg is the gluino mass, iriq is the squark mass and is the d'^ — h'^ mixing an- 

gle. The effective baryon number violating u'^d'^lf Yukawa coupling in the superpotential is 
parametrized here as f'{vR/Mp\) (see Eq. (2)). 

Let us first discuss the origin of the flavor mixing that changes a u'^d'^W operator to the 
required u'^d'^d^ operator. The dominant source for this in the present context turns out to 
be the mixing of h'^ with d'^. Such mixings occur in the left-right supersymmetric model since 
the right-handed quark mixings are physical above the scale vr. The renormalization group 
evolution of the soft SUSY breaking mass parameters between Mpi and vr will then induce 
mixings in the right-handed down squark sector proportional to the top-quark Yukawa 
coupling and the right-handed CKM mixings. This is analogous to the RGE evolution 
in the MSSM inducing squark mixing in the left-handed down squark sector proportional 
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to the left-handed CKM angles and the top-quark Yukawa coupling. We estimate this 
right-handed rf^ — b'^ mixing to be 

This estimate is obtained by integrating out the RGE between Mpi and vn assuming uni- 
versality of masses at Mpi |jT2|. Since above vji, both and b'^ are part of the same SU(2)ji 
multiplet, unlike in the MSSM, b'^ Yukawa coupling is of order one. In this momentum range, 
the top-quark Yukawa coupling reduces the mass of b'^. In going to the physical basis of 
the quarks, this effect will induce the squark mixing quoted in Eq. (5). For the numerical 
estimate we took mo = M1/2 and Aq = and vr ^ 2 x 10^^ GeV for illustration. 

There is a second source of flavor violation that induces d'^ — b'^ mixing in general SUSY 
models. That is the baryon number violating Yukawa couplings themselves. If we write in 
standard notation, the effective S-violating superpotential arising from Eq. (2) as W D 
{\'-ji^/2)u':d'jdl, the RGE evolution from Planck scale to the weak scale will induce d'^ — ¥ 
mixing proportional to ^'[jk- For example, if we keep only the couplings involving the u'^ 
quark, viz., A'/js, and \ii2^ we can estimate the induced (5)^^) by integrating the relevant 
RGE g to be 

Recalling that A" ~ vr/M^i, we see that while this source of flavor mixing may not be 
negligible, it would be typically smaller than the ones from the right-handed quark mixings 
ofEq. (5). 

A third source of flavor violation relevant foT N — N oscillations has been identified in 



Ref. [|14| involving the exchange of the Wino. Such diagrams will have an electroweak loop 
suppression and a chirality suppression necessary to convert the left-handed squark to the 
right-handed one. We find that this contribution to 5m^_^ has a suppression factor given 
approximately by [{a2/4:7r){mb/mg)]'^ ~ 1 x 10~^ (valid for small tan/5) which is about two 
orders of magnitude smaller in this class of models compared to the gluino exchange diagram 
of Fig. 1. 

One has to calculate the hadronic matrix element of the six quark operator in order to 



obtain the t^v-tv- This has been discussed in several places in the literature ||I5|. The calcu- 
lations of this "conversion" factor can be done using crude physical arguments, according to 
which one has to multiply the Gab=2 by |'?/'(0)|^ to obtain 5mj^_fi where is the baryonic 
wave function for three quarks inside a nucleon. On dimensional grounds, one can deduce 
that |^/'(0)|'^ ^ ^QCD' which implies that 5m^_f; ~ 10^^G'ab=2 GeV. More detailed bag 



model calculations have been carried out. Rao and Shrock in Ref. quote this conversion 
factor to be 2.5 x 10~^(jab=2- We shall use this number for our numerical illustrations. 
Combining this matrix element with Eq. (4)- (5) we obtain 



7 X lO^sec. /2 X 10^^ GeV^ / Mg 



^^-^^ /'2 { VR ) UoO Gevj UoO Gevj " 

We can rewrite Eq. (7) in a form that makes the connection with the neutrino mass 
more transparent. The mass of Vr can be expressed through the seesaw formula from Eq. 
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(2) as m^^ = {m^pYM-pi/ {fv"^ where m^p denotes the Dirac mass of z/^. Ehminating the 
high scale Vr from this, we have from Eq. (7), 



. 2.8 X lO^ee. J J J J J . (8, 

Since the value of m^^ can be determined from the atmospheric neutrino data under 
certain assumptions, we conclude that within the seesaw framework, measurement of N — N 
oscillation will be a measure of the Dirac mass of the tau neutrino. This can then be used 
as a way to discriminate between models of neutrino masses. 

To see the specific prediction for N — N oscillations within the context of the class of 
models under consideration, we need to know the z/,- Dirac mass. We can estimate it from 
the following relations for the Dirac masses of the third generation quarks and leptons in 
the SUSY left-right model: 

m = {ht^i cos au + ^f,2 sin 
rrih = {ht^i cos ad + /it,2 sin ad)Vd 
m^p = {h^- i cos + hr^2 sin a„)f « 
"rrir = {hrA cos ad + K,2 sin ad)vd (9) 

Here au,d are the Higgs mixing parameters obtained from Eq. (1) (Eg: tana„ = A1/A2) 
and are the VEVs of the MSSM doublets. From Eq. (9) it follows that in the limit 
oihti 3> ht2 and h^i ^ h^2^ "we get rriyD ~ m^ — . Since at such high scales ~ m,-, 
this predicts nij^p ~ m^. In fact, we find that unless the two terms in Eq. (9) for m^p are 
precisely canceled, the Dirac mass of z/^ will be approximately equal to rrit. 

Using m^p ~ mt and / ~ 1, /' ~ 10^^, we get a value for the N — N oscillation 
time which is tantalizingly close to the present experimental lower limit [§. For values of 
TTiyD 10 times smaller than m^, and taking the supersymmetric particle masses as large as 
1 TeV, we see that r^_jv is less than 9 x 10^ seconds,^ which is in the range accessible 
to a recently proposed experiment [|l^. It would thus appear that a search for neutron- 



antineutron oscillation will provide an enormously useful window into neutrino mass models 
and as such a powerful constraint on the nature of new physics beyond the Standard Model. 

The prediction for N—N oscillations can be sharpened if we make use of flavor symmetries 
to determine the coefficients / and /' in Eq. (8). We illustrate this with a specific choice 
of fiavor symmetry taken to be SU{2)h x U{1)h- The first two families of fermions 
form doublets of SU {2)h and have a [/(l)^ charge of +1 while the third family fermions are 
singlets under both groups. This fiavor symmetry is broken by a pair of doublets (f){l) +</){— 1) 
and singlets + xl"!)- Allowing for effective operators suppressed by a scale M larger 
than the VEVs of these fields provides a natural explanation of the fermion mass and mixing 



^We have not included the QCD evolution factor from the SUSY scale of few hundred GeV to the 
GeV scale. Based on the QCD factor of 1.33 for the three quark proton decay operator we 
estimate that the corresponding factor for the N — N case should be about 2, which would reduce 
the estimate of t^_7v by a factor of 2. 
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angle hierarchy. If we choose (0) = (^(pj = e^M and (x) = (x) = ^x^^ ^ resonable fit to all 
quark and lepton masses is obtained, including neutrinos, for e,^ ~ 1/7 and ~ 1/20 and all 



dimensionless couplings being order one |18]. In this model, we can estimate the couplings 



/ and /' from the horizontal quantum numbers. They are / ~ and /' ~ e^sin6'c, so 
thatQ ///'^ ~ 8 X 10"^ This estimate leads to t^v-tv — 2 x 10^ sec from Eq. (8). Allowing for 
uncertainties of order 1 in this estimate, we expect that Tj^^n ^^t to exceed about 10^° sec. 
Before we conclude a few comments are in order: 

1. The model becomes unacceptable as soon as SU{3)c x U{1)b-l is embedded into a 
higher symmetry such as SU{4:)c or SO(IO) group because in that case, the AB = 1 
operator described in Eq. (2) is accompanied by other R-parity violating operators 
coming from the same higher dimensional operator Oi^2 due to the higher symmetry. 
Together, they would lead to unacceptable proton decay rate. Thus observation of 
N — N oscillation would be a signal of an explicit SU{3)c x U{1)b-l symmetry all 
the way upto the Planck (or string) scale. 

2. Baryogenesis has to proceed through a weak scale scenario since the AB = 1 inter- 
actions in the model are in equilibrium down to the TeV scale and will wash out any 
primordial baryon or lepton asymmetry. We note that the baryon number violating 
interactions contained in O2 themselves can potentially be the source of weak scale 



baryogenesis ||19 



The lightest neutralino in this model is unstable and will decay via ~^ QQQ modes 
due to the presence of the effective AB = 1 operator This prediction is directly 
testable at colliders. An alternative candidate for dark matter must be sought. 



IV. CONCLUSIONS 

In conclusion, we have found that in a large class of seesaw models for neutrino masses, 
despite the high scale of sessaw dictated by the current neutrino oscillation data, neutron- 
antineutron oscillation is in the observable range. In fact, unless the Dirac masses of neu- 
trinos are far below those deduced under simple and reasonable assumptions, we predict 
an upper bound on the neutron-antineutron oscillation time in the range of 10^ — 10^° sec. 
This is very close to the present experimental lower limit on N — N oscillations. In the 
most conservative theoretical scenario, the measurement of N — N oscillation time would 
be a measure of the Dirac mass for the tau neutrino, given the values of squark masses. 
This in itself would be an extremely interesting result, since it would discriminate among 
theoretical models of neutrino masses. This is apart from the fundamental importance that 
any observation of baryon number violation will carry. We therefore strongly urge a new 
experimental search for neutron-antineutron oscillation. 



^Normally, the parameter / could have been of order one but in the horizontal model or Ref . ||T^ , 
due to large — i^r mixing, it is the flavor entry that dominates the atmospheric neutrino mass 
difi'erence and hence the horizontal suppression factor e^. The factor is due the fact that the 
operator must be invariant under U{1)h- 
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FIG. 1. Tree level gluino-squark diagram ior N — N oscillations. 
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